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Using a quantum mechanical model, the exact energy eigenstates for two-particle two- 
channel scattering are studied in a cubic box with periodic boundary conditions. A 
relation between the exact energy eigenvalue in the box and the two-channel ^-matrix 
elements in the continuum is obtained. This result can be viewed as a generalization of 
the well-known Liischer's formula which establishes a similar relation in elastic scattering. 
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1. Introduction 

Scattering experiments play an important role in the study of interactions among 
particles. In the case of low-energyhadron-hadron scattering, experimental results 
on the phase shifts are available l^. On the theoretical side, low-energy hadron- 
hadron scattering can be studied with a non-perturbative method, like Lattice QCD. 
In a typical lattice calculation, energy eigenvalues of two-particle states with definite 
symmetry are obtained. Therefore, it is important to relate the energy eigenvalues 
which are available through lattice calculations to the scattering phases which are 
obtained in the scattering experiments. This was accomplished in a series of papers 
by Liischer El for a cubic box topology. This formula, now known as Liischer's for- 
mula, has been utilized in a number of applications, e.g. linear sigma model in the 
broken phase^, and also in quenched El and unquenched^l QCD. 

For hadron scattering at low energies, elastic scattering dominates. However, 
when the energy of the scattering process exceeds some threshold, inelastic scat- 
tering starts to contribute. In the case of pion-pion scattering, for example, the 
scattering process is elastic below the four pion and the two kaon threshold. Al- 
though the four pion threshold is in fact below the two kaon threshold, four pion 
final states will not contribute significantly due to its weak chiral coupling to the 
two pion initial state. Experimental investigations also supports this argument. It is 
then interesting to study the relation between the multi-channel two-particle states 
and the scattering phases in general. In this work, we report a relation between 
the energy of a two particle state in a finite cubic box and the scattering matrix 
parameters. Further details are provided in Ref. 1^1. 



2 Chuan Liu, Xu Feng, Song He 



2. The Quantum Mechanical Model of Two Channel Scattering 

We study a non-relativistic quantum mechanical model of two-channel scattering 
with the following Hamiltonian: 



' ET-^V'j^W{r)V,{r)^ ■ 



with the potential vanishes for r > R. Et > designates a positive threshold energy 
of the second channel. 

Energy eigenstates of the Hamiltonian |^ with energy E can be decomposed into 
spherical harmonics: 5'i(r) = X); m where the radial wave-functions 

4'i;im{i') with i = 1, 2 satisfy the usual radial Schrodinger equation. Concerning the 
coupled radial differential equations, the following statement can be proven^: // the 
matrix valued potential V{r) is such that every matrix element ofr'^V{r) is analytic 
around r — and that limr^o r'^V{r) ~ 0, then the coupled radial differential equa- 
tions has two finite, linearly independent solutions near r — 0.' it[.*j(r), with i = 1,2, 

j — 1,2 such that: u|.'j(r) ~ r''Sij. 

At large r where the potential V{r) vanishes, the wave function of the scattering 
state can be chosen to have particular forms: 

s'''i-- + /ii(ki -f)^ 



This wave function has the property that in the remote past, it becomes an incident 
plane wave in the first channel with definite wave vector ki. It is an eigenstate of 
the full Hamiltonian with energy: E = kJ/(2TOi). Similarly, if the energy E > Et, 
one can also build another eigenstate of the Hamiltonian ^'(^^(r), which in past 
becomes an incident wave in the second channel with energy E — Et + Ic^/ {2m2). 

In partial wave analysis, one decomposes the coefficients: fij (scattering appli- 
tudes) into spherical harmonics. The scattering eigenstate in Eq. ||2Jl can aslo be 
decomposed accordingly with the radial wave-functions: 

^ 1 s[ie*''^^ + (— y+ig^^'^i'" 



2 q\')„i. 
7*^21 ^ 

and a similar equation for 5'(^^(r) and wf \r). It is obvious that the two radial 
wave functions w^^^ (r) and (r) thus defined form a linearly independent basis. 
They are hence linear superpositions of the general solutions: up^(r) and uf\r). 
The converse is of course also true. 

An important physical property is that the matrix elements which enter the 
expansion, namely ^^j^ form a 2 x 2 unitary matrix. In practice, this two-channel 
5- matrix is usually parameterized as: 
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where the parameters: 6{, 62 and rji are all real functions of the energy E. 



3. Energy Eigenfunctions on a Torus 

When enclosed in a cubic, periodic box with finite extension L, the Schrodinger 
equation of the system takes a similar form except that the potential is periodically 
extended and the eigenfunction has to satisfy the periodic boundary condition: [i?o + 
VL{r)]ip{r) — Eip{r) , tlj{r + Ln) — 'ip{r) . where Hq is the free Hamiltonian and the 
periodically extended potential is: VL(r) = J^n^i^ + L^)- The eigenvalue of the 
Hamiltonian now becomes discrete with smooth eigenfunctions. It is also convenient 
to partition the whole space into two regions. In the inner region, every point satisfies 
the condition: |r| < R, mod{L). In the outer region: f2 = {r| : |r| > R ,mod{L)}, 
where the interaction potential VL(r) = and the Schrodinger equation of the 
system reduces to two decoupled Helmholtz equations: (V^ + kf)ipi{r) = 0, i — 1,2. 
It it easy to see that: '^{r;E) = J^im I^Li ^im (^) ^m(n) solves the 

Schrodinger equation in the inner region for |r| < i? with being non- vanishing 
coefficients. In the outer region fi, the solution must be linear superposition of the 
singular periodic solutions of Helmholtz equation"^, we thus obtain a set of homo- 
geneous linear equations. In order to have non-trivial solutions for the coefficients, 
the corresponding matrix has to be singular. This condition then gives: 



(5) 



where the unitary matrices [/'■'■' are specific functions of the energy 1^. 

In most lattice calculations, the symmetry sector that is easiest to investigate 
is the invariant sector: . We therefore focus on this particular symmetry sector. 
In the first order approximation, if we neglect the mixing between the s-wave and 
(7-wave, we have for the sector: 

cos(Ai + A2-S°- 5^) = 7^0 cos(Ai - A2 - 5? + . (6) 

where we have also used the special parametrization Q for the s-wave S'-matrix 
elements. The quantities appearing above, namely S^, 62, rjQ, Ai and A2, are all 
functions of the energy: E = k'l/{2mi) — Et + A:|/(2m2). Eq. © is the simphfied 
formula for the s-wave S'-matrix elements when contaminations from higher angu- 
lar momentum (mainly from I = 4) are neglected. This relation is helpful since it 
provides a constraint on the four physical quantities. If the mixing of the o-wave 
is taken into account, the corresponding formula becomes more complicated"^. Al- 
though we worked out the formulae in a cubic box, similar relations can also be 
obtained for general rectangular box following the strategies outlined in Ref. ^3. 

Finally, let us speculate about possible exten sion to the case of massive field 
theory. As in the case of single channel scattering one can work in the Bethe- 
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Salpeter formalism for the two particle states and the results reported here can 
possibly be generalized to field theory. Details for this are now under investigation. 

4. Conclusions 

In this work, we have studied two-particle two-channel scattering states in a cubic 
box with periodic boundary conditions. Assuming that energy eigenstates are only 
two-particle states, the relation of the exact energy eigenvalues in the box and the 
physical parameters in the coupled channel S'-matrix elements in the continuum is 
found. This formula can be viewed as a generalization of the well-known Liischer's 
formula to the coupled channel situation (inelastic scattering). In particular, we 
show that the two-channel 5-matrix elements in the s-wave are related to the en- 
ergy of the two-particle system by a simple identity, if contaminations from higher 
angular momentum sectors are neglected. This relation is non-perturbative in na- 
ture and it will help us to establish connections between the S'-matrix parameters 
in the multi-channel scattering with the energy eigenvalues which are in principle 
accessible in lattice calculations. 
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